Manifestations of quantum holonomy in interferometry 
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Abelian and non-Abelian geometric phases, known as quantum holonomies, have attracted consid- 
erable attention in the past. Here, we show that it is possible to associate nonequivalent holonomies 
to discrete sequences of subspaces in a Hilbert space. We consider two such holonomies that arise 
naturally in interferometer settings. For sequences approximating smooth paths in the base (Grass- 
mann) manifold, these holonomies both approach the standard holonomy. In the one-dimensional 
case the two types of holonomies are Abelian and coincide with Pancharatnam's geometric phase 
factor. The theory is illustrated with a model example of projective measurements involving angular 
momentum coherent states. 

PACS numbers: 03.65.Vf 



I. INTRODUCTION 

The Abelian geometric phase in the sense of Berry 
[l| and Pancharatnam or non-Abelian holonomies 
in the sense of Wilczek and Zee kl are associated with 
curves in a Grassmann manifold i.e., the collection 
of all subspaces of a given dimension in a Hilbert space. 
Such curves may be realized in adiabatic evolution of 
a system dependent on external parameters [l], 0] or 
through a sequence of projective filtering measurements 
of observables [1, Q. In these contexts, non-Abelian 
holonomies arise in cases where the parameter depen- 
dent Hamiltonian is degenerate and where the measured 
observables have degenerate eigenvalues. The former sce- 
nario has attracted considerable attention in the liter- 
ature @, 0, !, US HO, El O and has recently been 
shown to be of relevance to robust quantum computation 
[3 OS El EE El, 03 HI HI HI- While the latter ap- 
proach to non-Abelian holonomies has been discussed in 
the limit of dense sequences of projection measurements 
in Ref. @ , a detailed analysis of the genuinely discrete 
non-Abelian setting, analogous to Pancharatnam's orig- 
inal discussion [2J of the Abelian geometric phase in the 
context of interference of light waves transmitted by a 
filtering analyzer, seems still lacking. 

In this paper, we examine quantum holonomy in 
the discrete setting, and thus complement the study of 
holonomies in the continuous setting pursued in Ref. [231 ] . 
We show that the discrete setting is "rich" in the sense 
that it admits more than one reasonable type of holon- 
omy. We demonstrate two distinct holonomies that arise 
naturally in this context We shall call these discrete 
holonomies 'direct' and 'iterative'. Although they are 
nonequivalent, the two types of holonomies nevertheless 
approach, in the limit of dense sequences, the Wilczek- 



Zee holonomy Q for closed paths, as well as its general- 
ization 23] for open paths, which appears to suggest that 
the extra richness of the discrete setting disappears in the 
continuous limit. Furthermore, in order to ensure that 
the direct and iterative holonomies are reasonable, we 
formulate them in terms of interferometric procedures, 
thus making them meaningful in an operational sense. 

The outline of this paper is as follows. In the next 
section, we introduce the concepts of direct and itera- 
tive holonomies in the Abelian case followed by their 
non-Abelian generalizations. We show how the two 
holonomies can be associated with the internal degrees 
of freedom (e.g., spin) of a particle in an ordinary two- 
path interferometer. Section IlIII contains an analysis of 
the case where one or several of the adjacent subspaces 
partially overlap, leading to the concepts of partial direct 
and iterative holonomies. An example involving sequen- 
tial selections of angular momentum coherent states is 
given in Sec. HVl The paper ends with the conclusions. 



II. HOLONOMY IN INTERFEROMETRY 

Relative phases can be measured in interferometry as 
shifts in interference oscillations caused by local manip- 
ulations of the internal states of the interfering particles. 
In its simplest form, this can be realized for a pure inter- 
nal input state ip that undergoes a unitary transforma- 
tion U in one of the interferometer arms. This results in 
an interference shift arg(-0|[/|V') and visibility |(V>|£/|'0 
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where the former is the Pancharatnam relative phase 

The above interferometer scenario can be used to de- 
velop two different holonomy concepts that are associated 
with the geometry of a sequence of points in a Grassmann 
manifold, i.e., the set of X-dimcnsional subspaces of an 
A^-dimensional Hilbert space. These concepts we shall 
call the direct and iterative holonomies of the sequence. 
The former type of holonomy is direct in the sense that 
the whole operator sequence representing the points in 
the Grassmannian is applied to the internal state in one 
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of the arms of a single interferometer. The latter type 
of holonomy is iterative in the sense that it is built up 
in several steps, where each step involves an interfer- 
ometer setup that depends on the preceding one. For 
one-dimensional (K = 1) subspaces, corresponding to 
sequences of pure states, the two holonomies are Abelian 
phase factors, while for higher dimensional subspaces 
(K > 1) they correspond to non-Abelian unitarities. In 
the following, we describe how the two types of quan- 
tum holonomies arise in interferometry in the Abelian 
and non-Abelian cases. 

A. Abelian case 

Let ipi, . . . , ip m be a sequence c of pure states with 
corresponding one-dimensional projectors • • •) 

\i/>m){i>m\- We assume that (i/j a+ i\ip a ) ^ 0, a = 
1, . . . ,m - 1, and (%l>i\ip m ) ^ 0. 

Let us first discuss the direct holonomy associated with 
the sequence c. Consider particles prepared in the state 

|*o> = |^>®^(|0) + |1)), (1) 

where \ip±) is the internal state, and |0) and |1) repre- 
sent the two interferometer arms. The internal state is 
exposed to the sequence of projection measurements cor- 
responding to c in the 0-arm, while a U(l) shift e iK is 
applied to the 1-arm. The filtering measurements corre- 
spond to the projection operators ir a — |?/; a )(?/'a| ( S ) |0)(0| + 
1 (g> |1)(1|, a = 1, . . . , m, where 1 is the identity opera- 
tor on the internal Hilbert space. A 50-50 beam-splitter 
yields the (unnormalized) output state 

!*(«)} = |(r[ c ] + e ^i)|Vi) |o) 

+i(r[c]- e < "i)|^>®|i>, (2) 

where F[c] = \ip m )(ipm\ ■ ■ ■ \ipi){ipi\- The shift of the in- 
terference oscillations in the 0-arm produced by varying 
k, is determined by the phase factor 

7 D = *[<V'i|r[ c ]|Vi)], (3) 

where = z/\z\ for any nonzero complex number z. 
The phase factor j£> is the direct holonomy of the se- 
quence c. 

The concept of iterative holonomy involves a sequence 
of interferometer experiments, each of which being de- 
pendent on the preceding one. Prepare the state 

|^ 1 ) = -^(|^)®|0) + |^)®|1)), (4) 

apply the U(l) phase shift e tK2 to the 0-arm, and let it 
pass a 50-50 beam-splitter to yield the output state 

l* 2 ' 1 ^)} = i( e ^ 2 > + |^>)®|0> 

+ \{e^ 2 )-\^))®\\). (5) 



The resulting intensity |||e lK2 |V'2) + l^i)!) m the 0-arm 
attains its maximum for e IK2 = e lK2 — $[(-02|'0i)]- Re- 
peat the procedure but with and e lK2 |-02) m Eq. 
(O replaced by e lK2 \ifj2) and e lK3 1-03), respectively. This 
yields 

+ |(e iK3 |-03> - e^|-0 2 >) ® |1> (6) 

and the corresponding interference maximum in the 
0-arm for e iK3 = e^ 3 = ^K^e^ 2 ^)] = 
<&[(V>3|V , 2)]3?[(V'2|V , i)]- Continuing in this way up to ip m 
and back to tpi results in the final phase shift 

e Sl = mHMM^m^m-l)} ■ ■ ■ $[(lh\4>l)]- ( 7 ) 

We define 77 = e iKl to be the iterative holonomy of the 
sequence c. 

Both 773 and 7/ are geometric in the sense that they 
arc unchanged under the gauge transformations \tp a ) — > 
e'*°|fj), a — 1, . . . , to, for arbitrary real- valued (3 a . Al- 
though operationally different, the direct and iterative 
holonomies 70 and 7/ are numerically equal. Indeed, we 
have 

7D = ^[(i>i\tp m )(i'm\lpm-l) ■ ■ ■ (V^l^l)] 

= ^[^il^^K^l^m-l)]---^^^!)], (8) 

which is 7/ according to Eq. |7|). In fact, 773 and 7/ are 
both equal to the Pancharatnam geometric phase factor 

BUS!. 

B. Non-Abelian case 

Consider a sequence C of discrete points pi,p2, ■ ■■ ,Pm 
in the Grassmann manifold, now with arbitrary subspace 
dimension K. There is a natural bijection between the 
Grassmann manifold and the collection of projectors of 
rank K. Thus, we may associate to C a sequence C 
of projectors Pi, ... , P m . We construct the intrinsically 
geometric quantity [23| 

r[c] = p m ...p l5 (9) 

which is the non-Abelian counterpart to T[c] in Eq. @. 
Physically, T[C] can be viewed as a sequence of incom- 
plete projective filtering measurements [f|. Let us in- 
troduce a frame T a = {|afe)}^L 1 for each subspace p a , 
a = 1, ...,m. The set of frames constitutes a Stiefel 
manifold, which is a fiber bundle [26[ with the Grass- 
mannian as base manifold and the set of i^-dimensional 
unitary matrices as fibers. We introduce the overlap ma- 
trix ni m 

{F a \Fb)kl = (a k \k), (10) 

which is used in Ref. p3| to define holonomy for a con- 
tinuous open path in the Grassmannian. The polar de- 
composition J {J- a \3~b)\U a ,b of the overlap matrix, where 
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[(^qIJ 7 ;,)! = yj {Ta\Tb){Tb\Ta), leads to the definition of 
relative phase U a ,b as 



U af> = \(T a \T h )\ {T a \T h ) 



(11) 



under the assumption that the inverse | (Ta\Tb) | exists. 
The existence of the inverse is guaranteed if l^al-^ft)! > 
0, in case of which we say the two subspaces p a and pb are 
overlapping (23[. For overlapping subspaces, the relative 
phase U a ,b is a unique unitary matrix. 
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FIG. 1: (Color online) Direct (upper panel) and iterative 
(lower panel) holonomy in the interferometer setting. In the 
direct scenario, the sequence Pi, ... , P m of filtering measure- 
ments is applied to the internal state in the 0-arm (upper 
path). In the 1-arm (lower path), the internal state is exposed 
to a unitary operation V. The intensity for each orthonormal 
basis vector of the initial subspace is measured in one of 
the output beams. Maximum of the total intensity, defined 
as the sum over all k, is obtained when V coincide with the 
direct holonomy of the sequence. In the iterative scenario, the 
internal states \a^) and |(a + are exposed to the unitary 
operations V a and V a +i, respectively, in the two arms. Maxi- 
mum of the total intensity (sum over k) is obtained by varying 
V a +i and keeping V a fixed. In this way, the unitary operators 
V2, ■ ■ ■ , V m , Vi are given in an iterative manner, yielding the 
iterative holonomy as the final unitary V\. 

We first consider the direct holonomy. A beam of par- 
ticles is prepared in an internal state represented by the 
vector |lfc) e T\ and divided by a 50-50 beam-splitter, 
yielding the state 
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|* fc ) = |l*)®-^(|0) + |l». 



(12) 



In the 0-arm the internal state is exposed to the se- 
quence C of projective filtering measurements, corre- 
sponding to the action of the projection operators Il a = 
Pa® |0> <0j + i® |1)(1|, a = 1, .. . ,m. A unitary V is ap- 
plied to the internal degrees of freedom in the other arm. 
The resulting state pass a 50-50 beam-splitter. The out- 
put intensity in the 0-arm reads 



ik = in 



(i*|rt[c]r[c]|i fc ; 



(13) 



where [V] H 

we have introduced the matrix product 

D = {Ti\T m ){T m \T m -x)...{?2\Ti) 
Summing over all k yields the total intensity 



lfc|V|lz) is a unitary K x K matrix and 

(14) 



K 



itot - 2_^ Ik - 4 

+ ^ReTr(y t £>) 



(15) 



Under the assumption that |_D| exists, the total inten- 
sity attains its maximum when 



D D. 



(16) 



The unitary matrix Up is the direct holonomy associated 
with the sequence C as measured in the interferometry 
setup shown in the upper panel of Fig. [T] 

Next, we consider the iterative holonomy, which, as in 
the Abelian case, involves the performance of a sequence 
of interferometry experiments. Suppose all adjacent sub- 
spaces of the extended sequence pi, . . ■ ,p m ,P\ are over- 
lapping. Prepare the state 



|^ 1 ) = ^(v 2 |2 fc )®|o) + |i fc )®|i; 



(17) 



where V2P2V2 = P%- A 50-50 beam-splitter yields the 
output intensity in the 0-arm as 



T 2,l _ 



life) +V 2 |2 fc )|| 

l + Re[{Fi\F 2 )V 2 ] 



/,/. 



(18) 



where [V2]., = (2fe | V2 12^) is a unitary K x K matrix. 
Summing over k yields the total intensity 

I 2 d = IX 1 = 1(k + BbTt[(^ 2 )V 2 ]), (19) 

fe=i 

which attains its maximum for V2 = V2 = t^2,i- In the 
next step, prepare 

2 > = 4| (^|3fe) ® |0) + V 2 \2 k ) O ]X>), (20) 
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where V 3 P 3 VJ = P 3 and (2 fc |U 2 |2 ; ) = [V 2 ] kV The two 
beams are made to interfere by a 50-50 beam-splitter. 
Adding the resulting output intensities yields 

^tot 2 = ^(l + ReTrfC/^K^ljPg)!^ 

x((7 3 , 2 t/2,i)V 3 ]), (21) 

which is maximal for V 3 = V 3 = U 3 ^U 2,1 ■ By contin- 
uing in this way up to P m and back to P\ , we obtain the 
final result 

V 1 = U I = C/i, m i7 m , m _i . . . U 2A . (22) 

The unitary matrix U 1 is the iterative holonomy asso- 
ciated with C. The interferometer setting giving rise to 
the iterative holonomy is illustrated in the lower panel of 

Fig. m 

Under the change of frames 

k K 
^«- > {X>k>[W r «] jfci ,} , o = l,...,m, (23) 
fc=i 

W a being unitary matrices, we have 

{F a +i\Fa) - wi +1 (T a+1 \T a )W a , 

U a +l,a — * Wl +1 U a+ha W a . (24) 

Such a change of frames can be seen as a gauge trans- 
formation, i.e., a motion along the fiber over each of the 
points pi , . . . , p m in the Grassmannian. From Eq. (|24p 

(7/ -> WlrZ/Wi, (25) 

i.e., the direct and iterative holonomies transform uni- 
tarily (gauge covariantly) under change of frames. 

The unitary matrices U d and U '/ are the non-Abelian 
generalizations of 7^ and 7/, respectively. However, 
while 7/3 = 7/, we have 17 jj 7^ C^z in general. There are 
situations, though, where the two approaches give the 
same result, e.g., for continuous paths in the Grassman- 
nian. This follows from the fact that for a smooth choice 
of 7^ = {\a k (s)}}f =v we have \(F s+5s \F s )\ = l + 0(6s 2 ), 
1 being the if-dimensional identity matrix. Thus, for 
s € [0, 1] we obtain 

D = (^ \^)(l + O(Ss 2 ))U hl - Ss ... 
x(l + 0(Ss 2 ))U Ssfi 
= {F \F 1 )Ui !l - Ss ...U Ss>0 + O(5s), (26) 

where the correction term is of order 0(Ss) since it con- 
tains <5s -1 terms. By using the assumption that |(7-o|7-°i)| 
is invertible and the fact that Ux^x—Bs . . . Us s ,o is guaran- 
teed to be unitary for sufficiently small 5s, we have \D\ = 

ICFol-Fi)! +0(Ss) and \D\~ l = \ (.Fo|.Fi)| _1 + 0(8s). It 
follows that 

U D = U 0A U lt i- Ss ...U Sst0 + O(6s) 

= Ui + 0(6s) (27) 



since |(7-o|7-"i)| (<7"b|<7"i) = Uq,i- Thus, in the Ss — > 
limit, we obtain 

U D = Ui = C/o.iPe-fo 1 A ^ ds (28) 

with [A(s)]fci = (cik{s)\ai(s)). In other words, in the con- 
tinuous path limit, the direct and iterative holonomies 
are equal to the Wilczek-Zee holonomy 3] for closed 
paths (for which t/0.1 = 1), as well as its generalization 
[23j for open paths. 

We finish this section by pointing out a relation be- 
tween the above iterative holonomy and the Uhlmann 
holonomy 1291 ] applie d to a special class of density oper- 
ators [lo, l3ll . [32I l33l. [H], [H, Hi| . This class consists of 
normalized rank K projectors, and we consider sequences 
jrPi, ■ ■ • , ■hPm, -jf-Pi 01 such density operators. If all the 
adjacent subspaces are overlapping, this is a sufficient 
condition for these density operators to constitute an ad- 
missible sequence 29] , for which the Uhlmann holonomy 
t/ u m reads 

C/uhl = Ux, m Um,m-l ■ ■ ■ ^2,1, (29) 

where J7 a +i, a , a = 1, • • • , m — 1, and Ux, m are the partial 
isometry parts of P a +iP a and P\P m , respectively [37l |. 
The overlap matrices can be written 

{T a +i\F a )ki = ^2((a+l) k \\P a+1 P a \\{a+l) n ) 

n 

x((a+l)„|tf" a +i 1 oh> ) (30) 

where |P a +iP a | is the positive part of P a +iP a . One can 
write (7-i|7-" m ) similarly. From Eq. (f30|) it follows that 
\U a +x,a]nl = {{a+l)n\Ua,+i, a \ai). By combining this with 
Eqs. and (gH}, and using that U a+1 . a = P a +iU a+1 , a , 
we find 

[Ui] M = (l k \U uhl \h), (31) 

for admissible sequences of density operators that are 
proportional to projectors [38l l39l| . 

III. PARTIAL HOLONOMIES 

If at least one pair of adjacent states in the extended 
sequence ipi,..., ipm , ipi are orthogonal, then the corre- 
sponding holonomies 7/j and 7/ are undefined. Simi- 
larly, U d and U 1 are undefined if any of the adjacent 
pairs of subspaces are orthogonal. However, the non- 
Abelian case includes partially defined holonomies, when 
the number of nonzero eigenvalues of the positive part 
of U \^ m U m ^ m -i . . . C/2,1 = / or D is greater than zero 
but less than the subspace dimension K. This occurs 
when at least one pair of adjacent subspaces is partially 
overlapping, which results in a nonunique unitary part of 
the overlap matrix. To remove this nonuniqueness, one 
may use the Moore- Penrose (MP) pseudo inverse [ioj . 
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denoted as G, to introduce a well-defined concept of rel- 
ative phase. Let T a and T b be two frames of two partially 
overlapping subspaces p a and p b . Then the MP pseudo 
inverse is obtained by inverting the nonzero eigenvalues 
of [(.Fal-Ft,)! in its spectral decomposition. We define 



invertible. In such a case, \D\ — \(J r \J r i)\ + O(Ss) 
and U i = Uq,iUi,i-Ss ■ ■ ■ Uss.o, where £/o,i is a partial 
isometry and Ui,i-$ s . . . Us s ,o is unitary. It follows that 



U D = Ui = C/oiPe^ 1 A ( s > 



(37) 



u a , b = \(f a \n)\ (fain) 



(32) 



as the relative phase between the two frames. The matrix 
U a ,b is a unique partial isometry. 

In Ref . [23| , the relative phase between frames of par- 
tially overlapping subspaces was used to introduce a con- 
cept of partial holonomy of continuous open paths in the 
Grassmannian. Here, we develop the corresponding con- 
cepts for the discrete sequence C. 

For the direct holonomy to be (totally) defined it is 
a necessary and sufficient condition that all the adja- 
cent subspaces (in the extended sequence) are overlap- 
ping [41] . Thus, if there is partially overlapping sub- 
spaces in the sequence, and if there is at least one nonzero 
eigenvalue of \D , then the MP pseudo inverse yields 



U 



D 



\D\ e D, 



(33) 



which we define to be the partial direct holonomy. 

In the iterative case, we again find that the holonomy 
becomes partial or undefined if at least one pair of adja- 
cent subspaces in the sequence C is partially overlapping. 
For such cases, I is not unitary since at least one of 
the matrices C/2,1, ■ ■ ■ , U\,m is a partial isometry. If |/| 
has at least one nonzero eigenvalue, we define the partial 
isometry part oil, i.e., 



U, 



in. 



(34) 



to be the partial iterative holonomy associated with C. 

Let us discuss how the partial holonomies behave under 
gauge transformations. From Eq. (|24p we obtain 



D 
I 



W\DW! 
W\lWi = 



\D\ -> Wjl-DlWi, 



(35) 



under the change of frames in Eq. (|23[) . Furthermore, for 
any matrix X and unitary matrices U and V we have 
(UXV) e = V^X e U^ (see, e.g., p. 434 in Ref. gj). 
Thus, 



\D\ 



W\\D\ e Wi, 
W\\l\ e Wx. 



(36) 



By combining Eqs. (|35|) and (|36|) . it follows that the di- 
rect and iterative holonomies transform unitarily (gauge 
covariantly) also in the partial case. 

We prove that the two partial holonomies in Eqs. (j3U|) 
and ([34"]) coincide with that of Ref. [23| in the continuous 
path limit. To do this, we consider the smooth choice 
T s = {|afe(s))}^ =1 and note that for sufficiently small 6s, 
the two holonomies become partial only if (Tq | T\ ) is not 



in the 5s — > limit, which is the partial holonomy of Ref. 



IV. ANGULAR MOMENTUM COHERENT 
STATES 

Consider a particle carrying an angular momentum j, 
j > 1. Let J na be the angular momentum component in 
the direction n Q characterized by the spherical polar an- 
gles 6 a ,4>a, i-e., n a = (sin9 a cos4> a ,sm9 a sin(j) ai cos9 a ). 
Let {|/i)}u=-j be the eigenbasis of J z . Consider a se- 
quence of filtering measurements of , a — 1, ...,m, 
each of which selects the maximal angular momen- 
tum projection quantum numbers /1 = ±j (angular 
momentum coherent states [42| ). The selection cor- 
responds to the two-dimensional projection operators 
Pn a = \j;n a ){j;n a \ + \ - j;n a ){-j;n a \, a = l,...,m, 
where | ± j; n a ) are eigenvectors of J na . The use of angu- 
lar momentum coherent states simplifies the subsequent 
calculation since \j; n a ) can be viewed as a product state 
of 2j copies of the spin-| state |i;n Q ), and | — j;n a ) 
similarly as 2j copies of | — |; n a ). 

Now, let (h = 1 from now on) 



T{6 a ,<Pa) = {e 



(38) 



For this choice of frames, the overlap matrix takes the 
form 



{Ho«Aa)\H6b,<t> b )) - {^ifls^by R(a!by 



(39) 



where 



R(a,b) 



S(a,b) 



a - O b 

cos 1 COS 



-I COS 



2 



2.; 



% + Ob 
-1 sin I — - — I sin 



2,7 



(40) 



We notice that {/\R(a,b)\ + {/}S(a,b)} = 1, i.e., the over- 
lap matrix cannot vanish for this system. 
If j is a half-odd integer, then 



(F(8*,<i>*)\H9b,cl>b)) 

= ^/\R{a,b)\> + \S{a,b)\ 2 U a , b , 



(41) 
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where U a .b is a unique unitary matrix. It follows that 
the direct and iterative holonomies are identical. 

When j is an integer, the overlap matrix in Eq. (|39p 
may have a nontrivial positive part. This implies that the 
two types of holonomies may be different. To illustrate 
this, consider the sequence of directions 111,112,113,114 
characterized by the polar angles (0o,</>o), (0i,0o), 
(0i,0i), (6o,<f>i), respectively. Assume that the first and 
third overlap matrices are degenerate. This happens for 
1 0i - O I = tt/2, which yields (F{0i,<M)\H0M) = 
2 1 ->U 2il and = 2 1 ^f7 4l 3. Here, 

(72,i = t/4.3 = |(l + Cx), where 1 and cr x are the 2x2 
identity and Pauli— X matrices, respectively. We further- 
more assume that 15(3,2)1 > |-R(3,2)| and |5(1,4)| > 



|i?(l,4)| for which a polar decomposition yields the uni- 
tary matrices 1/3,2 = e l:,XlCTz and U\^ — e^^ Xor7 % re- 
spectively, cr z being the Pauli-Z matrix. Here, Xk = 
2 arctan [ cos 0& tan (A0/2)] , k = 0, 1, where A</> = <j)\ — 
<f>o. We obtain the partial holonomies 



Ui 



qp l / e 



\q D \2\e^o 



qi 1 / e-^ xo 

M2 



e iV0 J 
e -UXo 



(42) 



where 



Vo 

<1D 
'II 



= — arctan 



1-sin 2 o sm 2 ^y sm{j X o) 
(l - sin 2 O sin 2 cos(jxo) + (-1)' (t 

A / \ .7 / 



sin 2 O sin 2 ^Y 



(l - sin 2 0i sin 2 ^Y cobOxi) + i^Y ^ <?i sin 2 , 
cos(jxi). 



(43) 



It follows that the direct and iterative holonomies differ 
unless qo and qi have the same sign and 770 = jxo- The 
latter happens only if sin 2 0o sin 2 (A</>/2) = 0. Note that 
U d {Ui) is undefined if qo — (qj = 0). 

V. CONCLUSIONS 

Corresponding to a sequence in the Grassmann man- 
ifold of if-dimensional subspaces in an iV-dimensional 
Hilbert space, we have defined two holonomies, which 
both are gauge covariant in the sense that they trans- 
form unitarily under the change of frames in the sub- 
spaces. Interferometer settings that give rise to the two 
holonomies have been delineated. In the non-Abelian 
case these two holonomies are generically nonequivalent. 
In the case of one-dimensional subspaces, however, both 



the holonomies reduce to the standard Pancharatnam 
phase. Moreover, we have shown that in the limit where 
the sequences form a continuous and smooth path in the 
Grassmann manifold, the two discrete holonomies coin- 
cide with the Wilczek-Zee holonomy Q in the case of 
closed paths, and its generalized noncyclic version [23[, in 
the open-path case. It is an interesting question whether 
there exist other discrete holonomies, distinct from the 
two considered here, that also converge to the standard 
holonomy in the limit of smooth curves, and if so, if those 
can be implemented intcrferometrically. 
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